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ABSTRACT 

Three cometary-shaped objects in the giant Hll region NGC 3603, originally 
found and identified as proto-planetary disks (ProPlyDs) by Brandner et al. 
(2000) using HST+VLT in the optical and near-infrared, have been detected with 
the Australia Telescope Compact Array {ATCA^) in the radio continuum at 3 and 
6 cm. All three ProPlyD-like objects are clearly resolved with an extent of a few 
arcseconds. The integrated 6 cm fiuxes are up to 1.3 times higher than the 3 cm 
fiuxes with spectral indices averaged over the whole clump between a = —0.1 and 
-0.5 {Si, oc I/"), indicating the likely presence of non-thermal emission in at least 
some of the sources. We present spectral index maps, and show that the sites of 
negative radio spectral indices are predominantly concentrated in the direction of 
the tails in at least two of the three ProPlyD-like nebulae while positive spectral 
indices are found in the region facing the ionizing star cluster. We propose that 
thermal bremsstrahlung and non-thermal synchrotron radiation are at work in 
all three ProPlyD-like sources. In at least one of the three objects optically 
thin non-thermal synchrotron emission appears to dominate when averaged over 
their whole spatial extent, while the spectrum of the second source shows a 
marginal indication of a non-thermal spectrum. The average spectrum of the 
third source is in agreement with thermal bremsstrahlung. All measured fiuxes 
are at least one order of magnitude higher than those predicted by Brandner et al. 
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(2000). Upper limits for mass loss rates due to photo-evaporation are calculated 
to be ~ 10~^ M0year~^ and for electron densities to be ~ 10^ cm~^. Due 
to the unexpectedly large radio luminosities of the ProPlyD-like features and 
because the radio emission is extended a (proto-) stellar origin of the non-thermal 
emission from a dust enshrouded star appears unlikely. Instead we propose that 
magnetized regions within the envelope of the ProPlyD-like nebulae exist. 

Subject headings: interstellar medium: Hll regions: individual (NGC 3603) — 
stars: pre-main-sequence — stars: early type, interferometry — stars: formation 



1. Introduction 

The giant Hll region NGC 3603, located at a distance of about 6 kpc, shows the densest 
concentration and largest collection of visible massive stars known in our Galaxy. Recent 
HST images of ~0'/2 angular resolution (Moffat et al. 1994) have shown that NGC 3603 

consists of three Wolf-Rayct (WR) stars and ~70 0-typc stars, with an estimated 40-50 of 
these stars located in the central ~30"x30" (1 pc x 1 pc) region. The three H-rich WR- 
stars of subtype WNL are also located within ~1" of each other in the central core and 
arc the brightest members of the cluster. They are believed to be massive main sequence 
stars which, like the bright stars in R136a, drive very strong winds (De Kotcr ct al. 1997; 
Crowther and Dessart 1998). NGC 3603 has also been shown to be a seat of active star 
formation outside the main cluster region (e.g. Brandner et al. (2000)). In particular, three 
distinct dense cometary nebulae were detected there using HST and VLT: Are these ultra 
compact H ll regions (UCHR) or ProPlyDs ? 

UCHRs are small (<0.1 pc) nebulae, internally photoionized by a deeply embedded, 
massive star. Even though the embedded stars are very luminous, they are invisible at 
optical wavelengths because of the surrounding dust; instead they show strong far- infrared 
emission. Their radio continuum emission is often associated with OH and H2O masers. 
The morphology of UCHRs in high resolution radio continuum images ranges from spherical, 
cometary, core-halo, shell, to irregular (Wood and Churchwell 1989). 

In contrast, a proto-planetary disk (ProPlyD) is a phenomenon describing a low-mass, 
young stellar object (YSO) with a circumstellar disk, embedded in a dense (neutral and 
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ionized) envelope that is being externally photo-evaporated by the ultraviolet radiation from 
one or more massive stars. These low-mass stars are not able to ionize their surroundings 
significantly. The disks in these systems are either observed directly (as in Orion where the 
disks are seen directly or in silhouette against the bright background nebula (Bally et al. 
2000)), or are inferred because in almost all of the ProPlyDs (young) low-mass stars are 
visible at optical or near-infrared wavelengths (see e.g. Stecklum et al. (1998)). At low 
angular resolution ProPlyDs look hke UCHRs, but at high resolution, optical and near- 
infrared images reveal their different nature. So most ProPlyDs would have been classified 
as UCHRs prior to those observations. In some cases the distinction between the two classes 
remains difficult. Given the fact that there are about four times more UCHRs than expected 
(Churchwell 1990) from star formation rates and the potential problems identifying host 
stars for UCHRs (e.g. some may not have enough FIR flux to account for an internal OB 
star), it may well be that a large fraction of catalogued UCHRs arc mis-identifled ProPlyDs. 
Interferomctric radio continuum observations of ProPlyDs arc needed to determine their 
structure, spectral indices and thus the nature of their emission, as well as mass loss rates 
and extinctions when compared to measured Ha fluxes. 

ProPlyDs were flrst identifled in the Orion Nebula^, where over 150 of them are known 
(O'Dell et al. 1993). Two ProPlyDs have been identifled in more distant nebulae: one 
in NGC 2024 by Stapelfeldt et al. (1997) and one, G5.97-1.17, in the Lagoon Nebula by 
Stecklum et al. (1998). The three ProPlyD-like objects, hereafter referred to PI, P2 and P3, 
in NGC 3603 (Brandner et al. 2000) are the biggest, youngest and most massive ones found 
so far. They are also the most distant known. These emission nebulae are clearly resolved in 
the HST/WFPC2 observations, and share the overall morphology of the ProPlyDs in Orion. 
All three nebulae are rim-brightened and tear-drop shaped with the tails pointing away from 
the central ionizing cluster. 

According to Brandner et al. (2000) the brightest object (PI), which has a projected 
distance of 1.3 pc from the cluster, has the spectral (excitation) characteristics of a UCHR. 
Optical spectra reveal the presence of an underlying, heavily reddened continuum source, 
which is also conflrmed by near-infrared VLT/ISAAC observations. The WFPC2 obser- 
vations show that only the outermost layer is ionized whereas the interior is neutral. The 



^There are a variety of different terms used in the literature for describing the knots of ionised gas in 
M42 (Orion): e.g. CKs (cometary knots), PIGs (partially ionized globules), EIDERS (external ionised 
(accretion) disks in the environs of radiation sources), and ProPlyDs (proto-planetary disks); for a summary 
see McCuUough et al. (1995). Another expression being used is EGGs (for evaporating gaseous globules) 
describing the objects found in M 16 (Eagle Nebula). The cometary knots found in the Helix Nebula are 
compact globules and very different from ProPlyDs as they contain no stars. 
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morphology of PI is described as a heart-shaped head with a coUimated structure in between, 
which can be understood as the superposition of two individual ProPlyDs. In contrast, P2 
and P3, located at projected distances of 2.2 pc and 2.0 pc from the stellar cluster, respec- 
tively, show approximately axisymmetric morphologies. No embedded disk or central star 
has been detected so far in any of these nebulae, preventing a clear identification as proto- 
planetary disks. The optical point source (see Fig. 2c) close to P3 is probably not physically 
linked to the nebula (Brandner et al. 2000). The ProPlyD-like structures in NGC 3603 are 
about two orders of magnitude fainter than typical UCHRs, but have a typical extent of 
9000 AU with tails extending to 21000 AU, much larger in size than the ProPlyDs in Orion. 

Recent 3.4 cm radio continuum and recombination line measurements of NGC 3603 
by De Pree et al. (1999), which were focussed on abundance measurements and the bright 
continuum emission from the ionized gas in this region, have an angular resolution of ~7" 
and a 5a sensitivity of 55 mJy. By obtaining high sensitivity and high angular resolution 
(~l"-2") ATCA observations we primarily aimed to study the radio emission of the winds 
from many of the early-type stars in the cluster, as well as to detect and resolve the ProPlyD- 
like objects and other gaseous regions in the cluster periphery. This paper will focus on the 
ProPlyD-like sources^ only, which have been detected and are shown to be clearly resolved 
with the ATCA. A subsequent paper containing a detailed study of the whole NGC 3603 
region based on our ATCA observations will follow. 



2. Observations and Data Reduction 

Radio continuum observations of NGC 3603 were made with the Australia Telecope 
Compact Array (ATCA) in the 6A, B, C and D configurations in five observing runs in 
February, April, June, September and November 2000 with a total of 60 hours assigned 
observing time. The observing frequencies were 4.8 GHz (6 cm) and 8.64 GHz (3 cm) with a 
bandwidth in each case of 128 MHz. Detailed observing parameters are given in Table 1. The 
three data sets were combined, reduced and analysed in MIRIAD using standard procedures. 
The full Stokes parameters were measured. The flux density scale was cahbrated using 
observations of the primary calibrator, 1934-638, assuming flux densities of 2.84 and 5.83 
Jy at 3 and 6 cm, respectively. 

A primary beam correction was carried out. To obtain high angular resolution and filter 



''Throughout this paper we use the term "ProPlyD-hke" for these cometary-shaped objects in NGC 3603, 
since the identification of these objects as true ProPlyDs is premature, owing to the lack of a clear detection 
of a central disk or star. 
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the extended emission to emphasize the small scale structure we used uniform weighting of 
the uv data while omitting the shortest baselines, up to 25 kA, respectively. We find that for 
a 25 kA cutoff all bright extended nebula emission is removed from the ProPlyD-like objects 
even at 6 cm. Fig. 1 shows an example of a combined 3 cm radio map where a uv cut of 
10 kA was chosen. A zoom on the three ProPlyD-like objects is shown in Fig. 2. For the 
6 cm maps a uv cut of 25 kA was used. The average r.m.s. in these combined radio maps was 
approximately 0.2 mJy at 3 cm and 0.4 mJy at 6 cm, but varies throughout the map. The 
r.m.s. levels near the ProPlyD-like objects, which lie in rather empty regions, were ~ 0.1 
mJy at 3 cm and ~ 0.2 mJy at 6 cm. The dynamic range in the 3 and 6 cm images is about 
1:100. Fitting a 2D-Gaussian to the dirty beam gives FWHM values of 0'/94 x 0'/77 to 1'.'06 
X 0'.'87 at 3 cm and 1'.'53 x 1'.'23 to 1'.'83 x l'/49 at 6 cm depending on the uv cut employed. 
In Fig. 2 we have restored the cleaned maps with a circular Gaussian beam of width 1" at 
3 cm and 2" at 6 cm. 

To create spectral index maps between 3 and 6 cm we convolved the 3 cm map with a 
Gaussian appropriate to achieve the resolution of the beam FWHM at 6 cm. The spectral 
index maps were then derived from the 6 cm and convolved 3 cm maps using a clip value of 
1.0 mJy beam~^, corresponding to approximately lOo" at 3 cm and 5a at 6 cm (statistical 
errors). Error maps for the spectral index maps were calculated using error propagation. — 
We note, however, that spectral indices of extended sources as obtained from interferometric 
data have to be regarded with caution. 



3. Results 

The three ProPlyD-like sources found with the ATCA at 3 and 6 cm are clearly resolved, 
showing a head-tail extent of ~4" (see Fig. 2). P3 shows the most pronounced head-tail 
structure with a 3 cm flux density ratio between head and tail of about 10:1. The tail is very 
well defined and at least 2" long, pointing away from the central star cluster. Unfortunately, 
P3 is rather faint in the low-sensitivity HST broad band image shown by Brandner et al. 
(2000); it is located outside the region of their high-sensitivity Ha image. 
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3.1. ProPlyD structure and fluxes 

The three ProPlyD-hke objects have a cometary or head-tail shape* at cm-radio wave- 
lengths, matching their appearance in optical images (see Fig. 2). 

Integrated radio flux densities and peak fluxes are derived by summing over the whole 
ProPlyD-like objects, and by fltting circular Gaussians to each object. Statistical flux un- 
certainties are estimated as map r.m.s. times number of beamsizes covered by the source. 
Dependent on the method of flux determination, beam size and uv cut employed to the maps 
(5, 10, 15, 20, 25 kA) the resulting flux densities vary. The standard deviation of this varia- 
tion is considered as systematic uncertainty. We found that the systematic uncertainties are 
always larger than the statistical error. In Table 2 we report the average fluxes and their 
uncertainties of all 3 nebulae, together with their average spectral indices derived from the 
integrated fluxes. 

Circular Gaussians^ are used to at least attempt disentangling the head and tail struc- 
tures in these cometary-shaped objects. The radio heads of the ProPlyD-hke objects are 
typically 2" -3" in diameter while the total head-tail extent is ~4". P3 shows the most 
pronounced head-tail structure of the three (most likely due to the viewing angle) and the 
smallest head (diameter '^2") among all three ProPlyD-like objects, with the largest peak 
flux at 3 cm while the integrated head flux is comparable. Its head-to-tail flux ratio is about 
10:1, whereas in the other two ProPlyD-like objects part of the tail may be confused with 
the head flux density. 

Both, PI and P2, appear inclined with respect to the plane of the sky, with PI inclined 
at a larger angle than P2. This can be seen in Fig. 2 where we have indicated the projected 
direction to the ionizing source. Fig. 3 shows the residual radio continuum image of PI 
after removing a single circular Gaussian component. The residual map of PI shows two 
excess emission sites, one corresponding to the tail, and a second one north of it. Recent Ha 
images indicate that PI is composed of two separate cometary-shaped objects (Brandner et 
al. 2000). We interpret the northern excess as the second head. A two component Gaussian 
flt to PI at 3 cm gives peak fluxes of 4.6 mJy beam~^ and 1.4 mJy beam"-*^, and integrated 
flux densities of 9.9 mJy and 5.7 mJy for the two heads, with a separation of the two peaks 
of 1'.'3. A two component model could not be applied to the lower resolution 6-cm data. 

The tails of all three ProPlyD-like objects are directed away from the star cluster. 



*The silhouette of ProPlyDs has variously been referred to as cometary, tadpole, tear-drop or simply 

head-tail shaped. 

^We find that elliptical Gaussians instead of circular ones would fit too much tail flux to the head. 
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approximately coincident with the tails in Ha within the positional uncertainty of the ra- 
dio/optical pointings (see Fig. 2; note the shift of 0'/5 in declination and right ascension 
of the HST-image which we have employed in these figures). The uncertainty of the radio 
position is ^C/.'IS while HST has a pointing uncertainty of < 1". 

The present data were obtained in all four Stokes parameters. No polarization signal 
above 3a was detected. 

Possible flux variability on timescales > days expected from e.g. flaring pre-main se- 
quence (PMS) stars with non-thermal radio spectra, cannot be tested with the present data. 

3.2. Spectral indices of the ProPlyDs 

The radio continuum fluxes (see Tabic 2) are a factor of 9 to 21 larger than those pre- 
dicted by Brandner et al. (2000) assuming optically-thin thermal bremsstrahlung as the radio 
emission mechanism (McCuUough et al. 1995), which scales with the (extinction corrected) 
Ha flux. The neutral part of the envelope of the ProPlyD-hke object may attenuate Ha 
photons from the far side of the object depending on the viewing angle. However, McCul- 
lough et al. (1995) have estimated this effect to be less than 25% of the total Ha light for an 
isotropic ensemble of ProPlyDs. Furthermore, the reasonable coincidence of the sizes of the 
ProPlyD-like objects in the radio and optical band suggests that this effect is not important 
here. For PI and P3 the integrated fluxes at 6 cm are higher than that at 3 cm, whereas 
similar flux densities were predicted. The average spectral index of each object, shown in 
Table 2, is calculated on the basis of the integral fluxes at 3 cm and 6 cm in this table. 
Modeling the uv distribution of the 6 cm map to match the uv distribution at 3 cm gives 
similar results. 

The average spectral indices of PI and P3 are negative (P3 only marginally), between 
a = —0.5 ± 0.2 and —0.3 ± 0.2 {Si, oc i/°), indicating non-thermal emission, while the 
spectrum of P2 (a ^ —0.1 ± 0.2) is in agreement with the expectation from optically thin 
thermal bremsstrahlung. 

Hydrodynamic simulations predict spectral indices of a ~ near the ionization fronts, 
which can increase up to a ~ 0.6 when considering in addition a hydrodynamically coUimated 
jet (S. Richling, private communication). Negative spectral indices cannot be explained by 
this model. 

Because of the high resolution and sensitivity in the radio maps, we were able to generate 
spectral index maps for all three ProPlyD-like objects. Fig. 4 shows examples of the spectral 
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index maps. Here we used the radio maps shown in Fig. 2 and constrain ourselves to 
brightness densities greater than 1.0 mJybeam"^, which corresponds to roughly 5 cr at 6 cm 
and 10 (T at 3 cm (based on statistical error estimates). The radio emission zone within 
each ProPlyD-like object appears to be inhomogeneous. The spectra in at least two of the 
sources appear steeper in the tail than in the head, with a tendency of positive spectral 
indices located towards the region facing the ionizing star cluster. 



3.3. Extinction towards the ProPlyDs 

Melnick et al. (1989) have derived the reddening towards NGC 3603 from UBV pho- 
tometry, and found Ay ^ A — 5 mag for the cluster core. However, as can be seen in their 
data, the visual extinction is not uniform across the whole star forming region, but increases 
significantly with distance from the star cluster, mainly towards the North and the South. 
This is confirmed by recent VLT/ISAAC data on the basis of 4750 stars, showing that the 
extinction can be as high as Ay ~ 14 mag (Brandner, private communication). For ex- 
ample, star MTT 68, close to P3, has Ay ~ 6.2 mag, whereas star MTT 81 near PI has 
Ay ~ 4.7 mag, and star MTT 79 close to P2 has extinction Ay ~ 5.2 mag. Here we use 
R = Ay/Eb-y ^ 3.2 as used by Melnick et al. (1989) for NGC 3603. Thus, it appears 
likely that P2 and P3 suffer from stronger foreground extinction than the stars in the cluster 
center which have on average Ay < 4.5 mag. 

Using the ratio between radio flux density and Ha flux as given by McCuUough et al., 
(1995; their equation 5) and assuming that the radio flux density at 3 cm is dominated by 
thermal bremsstrahlung, we predict the Ha flux using our 3 cm radio continuum data. The 
predicted fluxes are then compared to the Ha fluxes measured by Brandner et al. (2000) 
and used to derive lower limits for the extinction, Ahq. (If (gyro-)synchrotron radiation 
dominates the radio emission, which comes from an emitting region as specified in Sect. 4.4, 
the derived extinction An^ would be higher, with the exact value depending on the magnetic 
field.) Table 3 shows the calculated fiuxes (assuming — 10^ K) and extinctions. Note the 
uncertainty in comparing different areas for which the Ha fiuxes and the radio continuum 
fluxes have been obtained. Brandner et al. (2000) give the Ha surface brightness of the 
ProPlyD-like object heads, measured for a O'.'S aperture radius, whereas we quote peak 
fluxes per beam where the beam size is ~l'.'Ox I'.'O. For Ahq- ~ 0.85 Ay, we flnd Ay ^ 6.0 
for PI and Ay ^ 6.7 for P2. We flnd thus signiflcant excess extinction for both ProPlyD- 
like objects, which may be intrinsic to the ProPlyD-like objects themselves, otherwise the 
extinction estimate derived from the stars is underestimated. 
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3.4. ProPlyD densities and mass-loss rates 

Estimates for the electron/ion density inside the heads can be obtained from our radio 
data if thermal bremsstrahlung dominates the radio emission. This is hkely to be the case 
for P2 and maybe P3, while for PI the derived values should be considered as upper limits. 
The bremsstrahlung intensity in the radio domain for a thermal, ionized H-gas of constant 
density in a sphere is (Mezger and Henderson 1967) 

^radio ~ 1.9/^kpc/thCc^e,4iVion,4i^GHzr4"°-''niJy, 

where -D^pc is the source distance in kpc, ^sec is the apparent radius of the source in arcsec, 
iVe,4 and A^ion,4 are the densities of the thermal electrons and ions, respectively, in 10^ cm~^, 
z/GHz is the observing frequency in GHz and T4 is the electron temperature of the ionized gas 
in 10^ K. /th is the ratio of thermally emitting to total volume of the source, and is introduced 
to account for the possibility of smaller emitting volumes than the apparent extent of the 
source. This is hkely because only the outermost layers in ProPlyDs are beheved to be 
ionized. We estimate the ProPlyD heads ~ 2-3" in diameter. For T4 — 1 and A^e ~ -^ion 
we estimate the electron density by comparing the expected radio flux with the observed 
3 cm ProPlyD head flux densities, and find =(9, 7 and 13)-10Vth^^^ cm'^ for PI, P2 
and P3, respectively. (The corresponding emission measures are EM = 5.9, 3.5 and 9.3 
xl0^j\^^^^ cm~® pc.) If the emitting volume is smaller than the total source volume, the 
electron density and emission measure increase accordingly. From spectral considerations 
the emission must be optically thin, i.e. free-free absorption may set in at frequencies lower 
than 4.8 GHz. Thus values of EM cannot exceed ~ several 10^ cm~^ pc, in order to avoid 
free-free absorption at frequencies above 4.8 GHz. This constrains the size of the thermally 
emitting volume of the 3 ProPlyD-hke objects to at least 4%, 2% and 9% of the total source 
volumes, and consequently puts limits on the thermal electron densities in the sources: 4.5, 
4.9 and 4.3 -10^ cm~^ for PI, P2 and P3, respectively. Note that these numbers should be 
considered as upper limits if the observed radio flux is dominantly of non-thermal origin. 
These numbers are in reasonable agreement with Ne > 10'^ cm~^ derived from [SII] line 
ratios by Brandner et al. (2000). For the Orion ProPlyDs densities of order 10^ — 10^ cm~^ 
have been found (Henney and O'Dell 1999). Note that our derived numbers suffer from 
uncertainties due to the unknown geometrical projection angle of the ProPlyD-like objects, 
and hence the true source diameters. 

Mass loss from the ProPlyD-like objects occurs through external heating of the enve- 
lope by the far ultraviolet (FUV) radiation field from hot stars, especially from the core 
of NGC 3603, and subsequent evaporation. Knowing the electron density, the mass-loss 
rate in a wind is given by M = ATrR^NeiMnnVw, with R the system radius, uih the ion 
mass and vw the evaporation velocity. Assuming a pure hydrogen gas, /i — 1.3, we find 



-10- 



M/v2Q ~ (5 — 10) X 10~^ M0year~^ for the heads of all three ProPlyD-like objects, where 
^20 is the wind velocity in units of 20 km s~^. Brandner et al. (2000) found evaporation 
flow velocities of order 10-25 km s~^. The derived upper limits for the mass-loss rates are of 
the same order as estimated by Brandner et al. (2000), and about two orders of magnitude 
higher than the typical values for the Orion ProPlyDs. 

A rough estimate of the mass reservoir of the NGC 3603 ProPlyD-likc objects can tlion 
be derived by multiplying the mass loss rate with an estimated evaporation time scale of 
~ 10^ years (Brandner et al. 2000). This gives a mass for the ProPlyD-like objects of order 
1-10 M0, 10-1000 times larger than for the Orion ProPlyDs. The radius of the putative 
disks in these massive objects can be estimated using the evaporation model of Johnstone et 
al. (1997) which considers the mass reservoir in a disk-like shape. We find an approximate 
disk radius of 2000 AU (f« 0.3") for all three NGC 3603 ProPlyD-hke objects: Brandner et 
al. (2000) found ~3400 AU from hydrodynamical simulations. 

4. Discussion 

4.1. Comparison with the ProPlyDs in the Orion and Lagoon Nebulae 

The ProPlyD-like objects in NGC 3603 are 20-30 times larger and much more spectacu- 
lar than those in Orion, which were detected at 2 and 20 cm with the VLA (see McCuUough 
et al. (1995) and references therein; Henney and O'Dell (1999)). 

Orion contains numerous dense clumps, originally classified as UCHRs, but now known 
to be ProPlyDs. Most Orion ProPlyDs have central disks, sometimes seen only in silhouette 
against the background light. On recent HST-images dozens of jets powered by young stars 
embedded in ProPlyDs have been found (Bally et al. 2000). All Orion ProPlyDs detected at 
radio wavelengths are non- variable thermal emitters with spectral indices between -0.1 and 
0.2 (Felli et al. 1993), and located within 0.04 pc of the dominant ionizing 0-star [9^ Ori C). 
The non-thermal Orion radio sources show strong variability, and are mostly associated with 
visible pre-main sequence (PMS) stars. Their non-thermal spectra are usually explained as 
due to stellar flaring activity. 

G5. 97-1. 17 {D = 1.8 kpc) was originally classified as an ultracompact Hll region, but 
is now thought to be a ProPlyD at ~ 5000 AU projected distance from the 07 star Herschel 
36 in the center of M8, the Lagoon Nebula (Stecklum et al. 1998). ESO near-infrared, 
HST optical and VLA radio continuum observations indicate that G 5. 97-1. 17 is a young 
star surrounded by a circumstellar disk that is likely being photo-evaporated by Her 36, 
similar to the ProPlyDs in Orion. The previous hypothesis suggested that G5. 97-1. 17 is 
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a UCHR intrinsically ionized by an embedded BO star. However, Stecklum et al. (1998) 
showed that G5. 97-1. 17 is predominantly externally ionized, and thus the spectral type of 
the embedded star should be later than B5. The Ha flux over (/.'G is consistent with the 
A = 2 cm appearance. Optical and near infrared (NIR) continuum images show the central 
star is displaced from the peak of the bow shock by 0'.'125. NIR photometry of G5. 97-1. 17 
revealed that its central star is extremely red, which cannot be explained by extinction laws 
that use spherical matter distributions. This supports the idea of the circumstellar matter 
surrounding the central star of G5. 97-1. 17 being arranged as a disk. The cm-radio spectrum 
of G5. 97-1. 17 appears flat, which is interpreted as optically thin free-free emission (Wood 
and Ghurchwell 1989; Doherty et al. 1994). No OH or H2O masers, which are often associated 
with UGHRs, are known in G5. 97-1. 17. 

In comparison, the NGG 3603 ProPlyD-hke objects appear rather pecuhar and spec- 
tacular. Not only are they the largest and most massive ones found so far. Exposed to an 
extremely strong FUV radiation fleld, they also suffer the largest mass losses, though their 
distances from the ionizing source are larger than for the Orion or Lagoon Nebula ProPlyDs. 
No embedded source nor disk has been found for any of the NGG 3603 ProPlyD-like objects. 
At cm- wavelengths they are so far the only ProPlyD-like objects showing non-thermal radio 
emission. 

The spatial distribution of the NGG 3603 ProPlyD-like objects is also peculiar. They are 
not only distributed apparently along a straight line (sec Sect. 4.3 for further discussion), but 
all three also lie clearly to the north of the ionized gas region and avoid the region north-east 
and south-west from the stellar cluster. This is in contrast to the Orion ProPlyDs which do 
not appear to be located in preferred regions of the nebula (e.g. Bally et al. (2000)). Since 
sequential star formation proceeding from north to south seems evident in NGG 3603 (e.g. 
De Pree et al. (1999)), this suggests that the ProPlyD-like structures in this region might 
have emerged just recently. 

4.2. External versus internal ionization 

A possible way to prove that the ProPlyD-like objects have indeed been ionized exter- 
nally by the star cluster is to show that the observed radio fluxes for the ProPlyD-like objects 
are consistent with the number of ionizing photons they receive from the cluster stars. For 
simplicity we assume thermal radio emission, and note that this leads to an upper limit for 
the number of ionizing photons needed to produce the observed radio flux density. In this 
case the brightness temperature is directly proportional to the number of ionizing photons. 
Because NGG 3603 contains a large number of stars that are extremely hot and luminous 
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they provide an ionizing flux that is several orders of magnitude larger than that in Orion 
or M 8 (sec Tabic 4). The expected brightness temperatures at the projected distances of 
1.3 pc, 2.2 pc and 2.0 pc of the ProPlyD-likc objects from the ionizing source due to a Lyman 
flux of lO^^photons s^^ (Brandner et al. 2000) are 200 K, 70 K and 90 K at 3 cm assuming 
T4 = 1. This is only slightly higher than our estimates derived from the ATCA 3 cm maps 
of 40 K, 30 K and 70 K for the three ProPlyD-like objects, respectively. We can also com- 
pare the number of ionizing photons required to deliver the observed flux densities at the 
location of the ProPlyD-like objects, and compare this number with the number of Lyman 
photons they receive from the star cluster. Noting that the ProPlyD-hke objects subtend 
a fraction of (6, 3 and 2)-10~^, of the total sohd angle seen by the cluster stars, they thus 
receive (6, 3 and 2) -10^^ photons s~^, respectively. On the other hand, the required number 
of ionizing photons at the location of each ProPlyD-like object to deliver the observed radio 
flux density is about 2 • lO''^ photons s~^ in all three cases. These comparisons show that 
the number of required Lyman photons is slightly lower than provided by the star cluster, 
and thus would even allow an additional dust attenuation and/or, more likely, longer linear 
distances between the ProPlyD-like objects and the ionizing source in the external ionizing 
scenario. Signiflcant ionizing power from within the ProPlyD-like objects appears therefore 
unlikely, which puts the spectral type of any central stars to later than ~B1, and hence mass 
to below ~ 10 Mq. 



4.3. Do the ProPlyD-like objects host disks ? 

The present radio data as well as the recently pubhshed HST/VLT-images do not provide 
any direct evidence that the observed ProPlyD-like objects in NGC 3603 host any (proto- 
) stellar objects or disks. To a limiting Xg-magnitude of 18, no circumstellar disk or central 
stellar object in any of the three ProPlyD-like objects has been detected by HST-fVLT 
observations. Instead the cometary-shaped nebulae may rather be massive emission regions 
photo-ionized by the radiation from the star cluster. Mellema et al. (1998) have presented 
model calculations for such clumps, called FLIERs or ANSEA which are commonly located 
along the major axis towards the ionizing source. Indeed, all three ProPlyD-like objects 
apparently lie approximately along one straight line, with the outflow source, Sher 25, located 
roughly midway between PI and P3. Although the bipolar outflow from Sher 25 appears 
not to point in the direction of any of the three ProPlyD-like sources, the ring-like structure 
around Sher 25, roughly perpendicular to its jet, has been shown to expand with a velocity 
of ~ 30 km s~^ (Brandner et al. 1997), and could in principle be related to the origin of the 
ProPlyD-like clumps. In addition, models such as the FLIER model which are not based 
on disk-like mass distributions typically give short evaporation times, inconsistent with the 
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distance of the ProPlyD-like objects from Sher 25 (W. Brandner, private communication) 

and the approximate date of the mass-loss event about 6630 years ago in Sher 25 (Brandner 
et al. 1997). Thus, the existence of such features at large distances from the cluster centre 
appears to favor models where most of the material is concentrated in compact structures 
such as disks. The actual presence of disk-like structures is however not confirmed. 

4.4. The possible origin of the non-thermal emission 

The significantly non-thermal spectrum {a <^ 0, S^y oc z/°) from at least one of the 
three NGC 3603 ProPlyD-like objects (PI; P3 seems to be marginally non-thermal) appears 
to be steepest {a < —0.5) in the tail, slightly flatter {a = —0.5 to 0.1) in the head, and 
have a tendency of positive spectral indices {a > 0) located towards the region facing the 
ionizing star cluster rather than away from it. This is not expected from non-thermal radio 
emission from a wind-wind collision region. The non-thermal radiation must originate rather 
from inside the ProPlyD-like objects, either from their putative central stars hidden behind 
thick material, the surrounding disks or the envelopes. A stellar origin analogous to active 
star coronal emission seems to be unlikely because of the extended nature of the observed 
source. Non-thermal radio spectra could be produced by a population of energetic particles 
emitting synchrotron, gyrosynchrotron or gyroresonance radiation in a magnetic field. For 
gyrosynchrotron radiation the emission is concentrated at harmonic numbers s = v/vb = 
10 — 100 (where vb ~ 2.8-10^^i?G GHz is the gyrofrequency, Bq is the magnetic field strength 
in Gauss), and thus magnetic fields strengths of order 5 20 — 200 G are required. These 
high field values cannot be present in extended sources, which leaves synchrotron radiation 
as the only plausible emission mechanism. 

Disk or envelope, whatever the ultimate source might be, the negative spectral indices 
suggest optically thin non-thermal processes as the dominant radiation mechanism in PI, and 
maybe P3. One could argue that the shock created by the evaporation fiow predicted in Pro- 
PlyD models could accelerate charged particles within the objects. Other possibihties may 
be infiow-outfiow activity causing shocks during the early stage of star formation, compres- 
sion and reconnection of magnetic fields in the collapsing envelope, or magnetic reconnection 
in star-disk interactions in Young Stellar Objects leading to production of energetic parti- 
cles. Magnetic fields are known to exist in Molecular Clouds (Crutcher 1999), and therefore 
it appears plausible that gas clumps still contain "fossil" fields. In a magnetized plasma, 
acceleration due to particle collisions and subsequent bremsstrahlung can often be negligible 
in comparison with acceleration due to gyration around the field lines. In place of free-free 
emission, there is then synchrotron emission from energetic particles. The negative radio 
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spectral indices in PI and P3 are conveniently explained by non-thermal particle spectra. 
Because thermal bremsstrahlung is proportional to N^TfT^^'^ and synchrotron emission is 
proportional to Nf.B^ {j3 > 0), the former dominates if the density is high enough, or if 
the temperature or field strength is low enough. Indeed the electron densities found in the 
NGC 3603 ProPlyD-likc objects are ~two orders of magnitude smaller than in the Orion 
ProPlyDs (see Table 4). The different energy dependencies of the two kinds of emission may 
lead to one dominating at low frequencies and the other at high frequencies. 

A lower limit for a possible magnetic field is given by the Razin-Tsytovich effect. In order 
to avoid a Razin-suppression of the radio spectrum above A = 6 cm the magnetic field cannot 
be smaller than ~ 10~^ G for an electron density of 10"^ cm~^ in the emission region. On 
the other hand, to avoid gyroresonance absorption above A = 6 cm the field strength cannot 
be higher than 10^ G. The equipartition magnetic field, the minimum possible value in fiare 
loops, for a,Te — 10^ K particle distribution with density 10^ cm~^, is ~ 10~^ G. Note that 
synchrotron photons of energy 4.8 GHz in a 10~^ G field must have been produced by ~ 10 
GeV electrons. Efficient particle acceleration may occur, e.g., at the ProPlyD surface, from 
fiow amplification of even very weak ambient magnetic fields followed by local reconnection 
events (e.g. Shore & Larosa (1999)). 

The fiux density from a homogeneous sphere with diameter / due to synchrotron radi- 
ation in a uniform field and a low energy cutoff -Emin.McV in units MeV of a E~"'' particle 
distribution with ag = 2 spectral index^'' is given by (e.g. Dulk 1985) 

Su,syn ~ 5 X 10^^iVe_„t_4S^^£^min,MeV^pc-C'i7p^cZ^Gi/^ " 

The non-thermal electrons might have been energized by e.g. shock acceleration or magnetic 
reconnection events, out of the pool of thermal particles, or they may be Galactic cosmic 
rays penetrating the starforming region. In the latter case N^^nt ~ 10^^ . . . lO^^cm""' and 
ae ~ 2.4, and thus B ^ 10~^ G (for -Emm.MeV = !)• For the former case, because energy 
is conserved, the production efficiency of the high energy particles is related to the prop- 
erties of the thermal background particles, and vice versa. Realistic non-linear diffusive 
particle acceleration scenarios typically yield A^e,nt/-^e ~ 10~^ — 10"^ (i.e. the fraction of 
total electrons which end up with high energies; e.g. Ellison et al. (2000)), while magnetic 
reconnection might be as efficient as Nc^nt/Ne ~ 0.01 for the present densities and temper- 
ature (Tandberg-Hanssen & Emslie 1988). Since the thermal electron density cannot 
be significantly smaller than ^ lO^cm"^ (derived from [Sll] line ratios; see Sect. 3. 4), and 
-^e,nt ^ -^e if thermal and non-thermal electrons have the same source, we find that the size 



Shock acceleration typically predicts a, E ^ power law, which we use as a reasonable assumption here. 
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of the actual magnetized volume emitting non-thermal radiation must be smaller than the 
appearance of the ProPlyD-like object on the sky for this case. If particle trapping is at 
work, iVg nt > is possible, and consequently the ratio of emitting to total volume of the 
ProPlyD-like object becomes even smaller. 

The present data, however, do not allow to quantitatively constrain the actual size of 
the source responsible for the non-thermal radiation. 

We note that the recent discovery of variable X-ray emission from some ProPlyDs in 
Orion by Schulz et al. (2000) support the idea of a non-thermal emission volume smaller 
than the appearance of the ProPlyD-like object on the sky. The size of the emitting region 
was estimated on the basis of variability arguments to order 1-10 AU. 

5. Summary and Conclusions 

The three massive (~1-10 Mq) ProPlyD-like nebulae in NGC 3603, which have recently 
been discovered by Brandner et al. (2000), have been clearly detected and resolved with the 
ATCA at 3 and 6 cm, with one of them likely to be composed of two cometary-shaped 
objects. Their flux densities are about 10-20 times higher than expected from the Ha 
measurements, and a non-thermal average spectrum can be associated with at least one of 
the three ProPlyD-like sources. This is the first time that non-thermal radio emission has 
been detected from ProPlyD-like sources. Our spectral index maps show that the emission 
region is rather inhomogeneous, with negative spectral indices in the tail and part of the head 
whereas positive spectral indices, indicating thermal free-free emission, tend to be detected 
from a small region facing the star cluster. 

We derive upper limits for the mass-loss rates of 10~^ M© year~^ and electron densities 
of 10^ cm~^. These are in reasonable agreement with estimates from recent HST-images. 

We show that synchrotron emission from a magnetized, relativistic, non-thermal particle 
population may explain the non-thermal spectral regions. Energetic electrons necessary for 
synchrotron radiation may be produced through shocks or by magnetic rcconnection out of 
the pool of thermal particles, or may be Galactic cosmic rays. 

A stellar origin of the observed flux densities appears unlikely considering that the 
sources are extended. Thus magnetic fields, which have been neglected in ProPlyD models 
so far, are possibly associated with the ProPlyD envelopes or disks, and appear to be crucial 
in understanding the physical processes in these ProPlyD-hke objects. 

Melnick et al. (1989) have shown that an inhomogeneous dust distribution in NGC 3603 
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causes increasing extinction with distance from the star cluster. In particular, the extinction 
at the location of the ProPlyD-like objects is estimated to lie between Ay ^ 5 — Q mag. 
Extinction estimates derived from the radio- to Hct-luminosity ratio are 1-2 magnitudes 
higher. 

The ultimate discovery of the so-far undetected disks, which are thought to be part of 
all ProPlyDs, would give important hints about the origin of the extremely high and non- 
thermal radio fluxes and the nature of the comet ary- shaped clumps seen in NGC 3603. This 
may be accomplished in future mm observations with an upgraded ATCA. 

We are grateful to Wolfgang Brandner for providing the HST image of NGC 3603 to over- 
lay on our radio images, for interesting discussions and carefully reading the manuscript. We 
thank Jessica Chapman, Jim Caswell, Simon Johnston, Tylor Bourke and Sergey Marchenko 
for reading the manuscript and many constructive comments which improved our paper sig- 
nificantly. AM acknowledges a postdoctoral bursary from the Quebec Government. AFJM 
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Fig. 1. — 3-cm continuum map with a uv cutoff at 10 kA of NGC 3603 overlaid onto the 
Ha and broad band HST-image from Brandner et al. (2000). Note that in this radio map 
the small scale structure is emphasized due to our chosen uv-cutoS in the data analysis. For 
smoothing purposes the map was restored with a 2"x2"beam size. The brightest, extended 
continuum regions correspond to the heads of the giant gaseous pillars, showing evidence for 
the interaction of ionizing radiation with cold molecular hydrogen clouds. The radio emission 
from IRS 9, a deeply enshrouded association of protostars at the foot of the South-Eastern 
pillar, will be discussed in a forthcoming paper (Miicke et al., in prep.). The three nebulae 
indicated are the ProPlyD-like sources; their tadpole shape is shown in Fig. 2. The contour 
levels at 3 cm are -0.8, 0.8, 2.0, 4.0, 6.0, 12.0, 24 and 48 mJybeam-^ The beam is shown 
inside the box at the bottom left corner. North is to the top and East to the left. 
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RA (J3000) RA (J3000> RA [J30D0) 



Fig. 2. — ATCA 3 cm (upper row; uv cutoff at 10 kA) and 6-cm (lower row; uv cutofT at 25 
kA) radio continuum contour maps of the three ProPlyD-hke sources, numbers 1-3 from left 
to right, overlaid onto the Hct image from Brandner et al. (2000) for PI and P2, and onto 
the HST broad band image for P3. The contour levels for the 3 cm maps are -0.45, 0.45, 
0.75, 1.05, 1.5, 2.25 and 4.5 mJybeam"^, and for the 6 cm maps are -0.9, 0.9, 1.5, 2.1, 3.0, 
4.5 and 9.0 mJybcam^^. The beam is shown in the bottom left corner of each map. The 
Ha images have been shifted 0'.'5 to the North and 0'.'5 to the East in better alignment with 
the radio images. The arrows indicate the directions projected onto the sky from which the 
stellar cluster winds come. The cross near P3 indicates the position of the Chandra X-ray 
source. 
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Fig. 3. — ATCA 3 cm residual image of PI after subtracting a single circular Gaussian 
overlaid onto its unsubtracted brightness contour map as shown in Fig. 2a. The contour 
levels correspond to -0.3, 0.3, 0.5, 0.7, 1.0, 2.0 and 5.0 mJybeam"^ PI appears to possess 
two heads. 



RA (J3000) 



RA (J2000) 



Fig. 4. — Examples of radio spectral index maps {a; colour-coded) with the corresponding 
uncertainty levels (Aa; contours) for PI, P2 and P3 (from left to right) derived from the 
radio maps shown in Fig. 2. The contour levels are Aa = 0.2, 0.3, 0.4, 0.5 and 0.6. The 3 
and 6 cm flux densities used to produce these maps have been clipped at 1.0 mJybeam"^. 
Note the large spectral index uncertainties at the rim of the ProPlyD-like objects due to the 
low flux densities there. 
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Table 1: Observing parameters. 



Date (local time) 
Configurations 
Time on source 
Pointing position 
RA (J2000) 
DEC (J2000) 
Total bandwidth 
No. of channels 
Frequencies 
Beam size 

average r.m.s. near 
ProPlyDs [mJybeam 
Flux calibrator 
Phase cahbrator 



-n 



February 9/10 April 8/9 June 19 Sept. 13/14 Nov. 8/9 
6A 6D 6B 6A 6C 

561 min. 514 min. 480 min. 513 min. 543 min. 

11*^15"^ 07' 
-61° 16' 00" 
128 MHz 
32 

8640 MHz, 4800 MHz 
~l"x 1"(3 cm) 
~2"x 2" (6 cm) 
-0.1 (3 cm) 
~0.2 (6 cm) 
1934-63 
1059-63 



Table 2. Position, fluxes and spectral indices of the ProPlyD-like objects. 



Name 
A 


PI 

3 cm 


6 cm 


P2 

3 cm 


6 cm 


P3 

3 cm 


6 cm 


Radio position^ 


11:15:13.042, 
-61:15:50.38 


11:15:13.058, 
-61:15:50.53 


11:15:16.557, 
-61:16:06.57 


11:15:16.579, 
-61:16:06.62 


11:15:01.336, 
-61:14:45.69 


11:15:01.330, 
-61:14:45.53 


Peak fliix 

(mjy beam~^) 


4.3 ±0.5 


11.5 ± 1.6 


3.2 ±0.4 


7.8 ± 1.1 


5.6 ±0.7 


10.9 ± 1.4 


Integrated 
flux^(mJy) 


13.9 ± 1.1 


18.1 ±1.2 


11.3 ±1.4 


12.1 ±0.7 


14.1 ±1.0 


16.5 ±0.7 


Spectral index 
of whole source^ 


-0.5 ±0.2 


-0.1 ±0.2 


-0.3 ±0.2 


JHa''(cm-2s-i) 


0.56 


0.56 


0.20 


0.20 






Predicted 
radio flux^(mJy) 


1.56 ±0.39 


1.65 ±0.41 


0.56 ±0.14 


0.59 ±0.15 







Note. — The position of P3 given by Brandner et al. (2000; their Table 1) contains a typographical error. 

^Positions correspond to the radio locations of the ProPlyD-like object heads. 

^The quoted 3 cm and 6 cm integral fluxes are derived from images with different resolution. 

^Flux density Sv oc v", a = spectral index. 

^Ha flux corrected for an assumed foreground extinction in Ha of Ana = 4 mag. 

^See Brandner et al. (2000): assuming an electron temperature for the Donised gas of 10^ K; the flux uncertainty 
is less than 25% due to attenuation by the neutral envelope of the ProPlyD-like object itself as estimated by 
McCullough et al. 1995 (see text). 
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Table 3: Fluxes and derived extinction 



ProPlyD 




1 


2 


3 


3 cm peak flux 


[mJybeam ^] 


4.3 


3.2 


5.6 


predicted Ha peak flux 




1.5 


1.1 


2.0 


derived A^.^ 


[mag] 


5.1 


5.7 
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Table 4: Comparison with other known ProPlyDs^^^ 



Orion Nebula Lagoon Nebula 

M42 M8, NGC 6523 NGC 3603 



distance 


[kpc] 

L -L J 


0.45 


1.8 


6.1 


central star(s) 




el Ori C 


Herschel 36 


cluster 


star type(s) 




07 


07.5V 


3 WNs+ ~70 stars 


-^bol 


[L©] 


~ 10^ 


~ 10^ 


2 xlO^ 


Lyman flux 


[photons s~^] 


~ 8 • 10^^ 


2 xlO^« 




extinction 


[mag] 


4-6 


~ 5 


4-6 


masers 




H2O, SiO 


none known 


OH, H2O, CH3OH 


number of ProPlyDs 




> 150 


1 = G5.97-1.7 


3 


head size 


[AU] 


45-355 


1080 


7200-10800 


distance to central stars 


[pc] 


0.01-0.15 


0.024 


1.3, 2.2, 2.0 


radio flux 


[mJy] 


0.3-1.4 


17 


10-12 


brightness temp. 


[K] 


~ 10^ - 10^ 


~ 1500 


30-90 


emission measure 


[pc cm~^] 


S-IO^ 


(0.3 - 2) • 10« 


(3 - 9) • lOVt^"^"^') 


electron density 


[cm~3] 


10^ - 10^ 


(4 - 20) • 10^ 


~ 10^ 


mass loss rate 


[Moyr-i] 


1.2 • 10-^ 


7 ■ 10-^ 


10-5 


bow shock distance 


[AU] 


~ 100 


540 


44400,17400, ••■ 


velocity 


[ kms~"^] 


~ 10- 15 


~ 10 


10-25 


evap. time scale 


[years] 


~ 10^ 


~ 10^^ 


10^ 


ProPlyD mass 


[Mo] 


0.01-0.1 


~ 0.1 


~ 1 - 10 


disk radius 


[AU] 


27-175 


160 


3400 


disk mass 


[Mo] 


~ (0.5 - 20) • 10-3 


~ 6 • 10-2 


not known 



(^^Dctailed information on the NGC 2024 ProPlyD is not available. 
(2) See Sect. 3.4. 

References. — Stecklum et al. 1998, Brandner et al. 2000, Henney & O'Dell 1999, McCuUough et al. 
1995 



